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m ■ Abstract 

o, 

We review the recent results in the physics of the magnetic monopoles in gluo- 
dynamics and a dual formulation of non-Abelian theories, relevant to the physics of 
the confinement. It occurs that the dual gluon is a U(l) gauge boson, despite of the 
fact that usual gluons are non-Abelian. The effective infrared Lagrangian for glu- 
odynamics is suggested which leads to the Casimir scaling of the string tension for 
quarks in various representations. We also show that the results of the calculations 
in lattice gauge theories confirm our theoretical predictions. 



1 Introduction 



Magnetic monopoles were introduced theoretically by Dirac about 70 years ago 0. By 
definition, the radial magnetic field of the monopole is similar to the electric field of a 
point-like charge: 



H = 7ZZ3 > (i; 



> Coul Q r 

4:7rr' c 

where Q m is the magnetic charge. However, the field ([I]) by itself is inconsistent with 
the Maxwell equation divH = 0. To ensure the conservation of the magnetic flux Dirac 
postulated, therefore, that the magnetic field is transported to the monopole from infinity 
through a string: 

H = u. Coul + H Str , H str = -QzS(x)8(y)e(z), (2) 

where the Dirac string is directed along the negative z axis and z is the unit vector. The 
corresponding monopole gauge potential has a simple form in the spherical coordinates: 
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A = A,= , ^ = ^<i±^>. (3) 

4f r sin 6 

Thus, monopoles cannot be introduced without a string attached to it and, at first 
sight, any similarity between the electric and magnetic charges is lost. The only way to 
nevertheless maintain the similarity between the two types of charges is to ensure that the 
Dirac string is in fact unobservable. Thus, Dirac derived a few constraints on the theory 
which should be satisfied to make the monopoles effectively particle-like, not string-like 
objects. Theoretically, the Dirac monopole and the Dirac strings are fascinating subjects, 
for review see, e.g., [|[ §]. Still, they look exotic objects with infinite energy and have 
never been observed. 

The saga of the monopoles took a new turn with formulation of the (dual) supercon- 
ductor model of confinement @. It is well known that the color states are not observable 
as free particles since the potential between two static quarks grows linearly with the 
distance r at large r : 

lim Van(r) = a ■ r 

1 — >oo 

Note that the potential V qq - is defined in a gauge invariant way in terms of the Wilson 
loop Wc (for a pedagogical discussion and further references see, e.g., Ref. ||). Namely, 

V qq = - lim i log < W Fo >, (4) 

1 — >oo ± 



where 



Wc = TrP exp \ i <j> A^dx^j (5) 



and is the non-Abelian gauge potential while T is a rectangular contour r x T. 

Searching for an analogy to this spectacular phenomenon of confinement, one can 
notice that if magnetic monopoles, as test particles would be placed into a superconductor 
the heavy monopole potential would grow linearly as well. Indeed the magnetic field 
cannot penetrate the superconductor and would stream into an Abrikosov vortex . The 
superconductor model of confinement (for review see, e.g., |7|, [9]) assumes that the color 
electric field of the heavy quarks is organized into a similar tube-like structure in the 
vacuum state of the gluodynamics. For the model to realize, one needs condensation of 
the magnetic monopoles in the QCD vacuum. 

Although the model looks very appealing and, in a way, no alternative to it has ever 
been found, it poses quite a few questions which seem very difficult to answer. First of 
all, we consider now a non-Abelian theory while the (dual) superconductor model of the 
confinement copies the electrodynamics, that is a U(l) theory. Moreover, there are no 
stable classical monopole solutions to the pure Yang-Mills equations with a finite energy || 
and the condensation of the monopoles is difficult to visualize. 

The scope of this review does not allow even to mention the milestones on the road to 
the resolution of the puzzles brought by the model {§. To make the story very short, the 
monopoles were copiously observed and the dual superconductor model of confinement 
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has been confirmed, for review see, e.g., f|. However, there was a price to be paid. Just 
to mention a few points: the lattice regularization and numerical methods seem crucial, 
as well as the choice of a particular gauge. 

The need for the lattice regularization in the ultraviolet might look most surprising 
since the motivation to introduce the monopoles was to explain the confinement, i.e. the 
basic feature of the physics in the infrared. However, there are good reasons for this. 
Indeed, the Dirac strings attached to the monopoles are infinitely thin in the contin- 
uum limit. Therefore, the Dirac constraints on the strings can be checked only within a 
particular ultraviolet regularization. Moreover, the gauge potentials associated with the 
monopoles are singular both at the string and at small r, see, e.g., Eq. (|3|). In other 
words the monopoles, also in non-Abelian theories represent topological defects and the 
potentials are to be regularized in the ultraviolet. 

Moreover, the use of the numerical simulations is not only due to the complexity of the 
vacuum state but, even more so, due to the lack of understanding of a single monopole 
configuration. Indeed, classically only the monopoles with the minimal magnetic charge 
(in the units of the Dirac quantization condition) are stable and have an infinite action |J . 
The "lattice monopoles" observed in the numerical simulations would correspond to un- 
stable classical filed configurations and cannot be, therefore, understood quasiclassically. 
In this respect the lattice monopoles in non-Abelian theories are very different from in- 
stantons in case of the same non-Abelian theories or Dirac monopoles in the compact 
17(1) 0. 

The use of the specific gauges also appears as a precondition for the success of the 
dual superconductor model. Partly, this can be understood theoretically. Indeed, it was 
realized long time ago (Tl[ that the monopoles can emerge as effective dynamical degrees 
of freedom only if one fixes gauge up to a remaining U(l) symmetry. However, the 
numerical simulations indicate strongly that various ways of fixing the remaining U(l) 
are not equivalent to each other at all as far as the phenomenology is concerned. The so 
called Maximal Abelian gauge turns to be most successful |L2], ||. 

Schematically, the idea of the monopole condensation was formulated to explain the 
basic feature of the QCD as we know it in the real world, that is confinement. The idea 
was brought to the form suitable for the lattice simulations and dramatically confirmed 
through such simulations. However, the monopoles observed on the lattice look quite 
specific just for the lattice formulation and we need now to address anew the continuum 
theory to incorporate the lessons learnt through the simulations. One can say that the 
pendulum is swinging now in the opposite direction: from the lattice formulation towards 
the continuum theory. Broadly speaking, this review is summarizing the recent develop- 
ments in this direction. While the review is based mostly on the original papers in Ref [[0| , 
it is worth emphasizing that there are other papers devoted to the same or related topics, 



see, e.g., [fL4|, jig , jig] . Thus, it appears quite a common trend and we believe that work in 
the direction of incorporating lessons from the lattice studies into the continuum theory 
would continue in the future. 

More specifically, we will concentrate on the question how to bridge the U(l) nature 
of the monopoles with the observation that they play a key role in the dynamics of non- 
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Abelian theories. This double-face nature of the monopoles is revealed most straightfor- 
wardly through study of the heavy monopole potential V mm (r), or the so called 't Hooft 
loop ]17| . The 't Hooft loop allows to introduce a pair of (infinitely heavy) external 



monopoles, like the Wilson loop ([5]) allows to introduce a pair of external (infinitely 
heavy) quarks. It is worth emphasizing already at this early point that the monopoles 
introduced via the 't Hooft loop have the minimal magnetic charge Q m = 1 in the units 
of the Dirac quantization condition while the lattice monopoles living in the vacuum have 
Q m = 2. In particular, the Dirac string corresponding to the Q m = 1 monopoles is not 
visible to the gluons but would be visible to quarks while the Dirac string attached to 
the Q m = 2 monopoles is not visible both to particles in the adjoint and fundamental 
representations of the color group. 



The 't Hooft loop operator is most easily formulated on the lattice |18| . The monopoles 



are understood then as the end-points of the Dirac strings which in turn are defined as 



piercing negative plaquettes (see Section |3]2| for a detailed discussion). In the language of 



the continuum theory it would be natural to define the 't Hooft loop in terms of the dual 
gluon field 



H c = exp y~J c B v dx vj 



(6) 



Then the heavy monopole potential is defined similarly to 



V mm = - lim -log< H To >. (7) 

1 — >oo ]_ 

The central point is of course, how to relate the "dual gluon" to the fields entering the 
standard Lagrangian of the non-Abelian gauge theory. And this is one of the central 
issues which we will try to clarify here following ]13| . 



Generally, a dual gauge boson is understood as a field interacting with the magnetic 
current j m : 



J int 



Qm,B " jV, 



One of our basic points is that the field B should be treated as fundamental in gluody- 
namics since the external monopoles introduced via the 't Hooft loop are point like in 
the continuum limit. A well known example of introduction of a dual gauge boson is the 
Zwanziger Lagrangian describing electrodynamics with both electric and magnetic point 
like charges fl!9| . Formally one introduces two vector fields, the "standard" photon 



and the dual photon B^, so that L int = Q e A ■ j e + Q m B ■ j m . However, the number of the 
degrees of freedom is not changed since there is a constraint that the field strength tensor 
constructed on the potential A coincides with the dual field strength tensor constructed 
on the potential B^. More precisely: 

mpF^B) = m^F^A) (9) 
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where m M is an arbitrary space-like vector. The choice of the vector m M is a kind of 
new gauge freedom. Physically, one can visualize m M as the vector along which all Dirac 
strings are directed. 

In case of the gluodynamics, the dual gluon is still a U(l) gauge boson despite of 



the fact the direct gluons A a are in an adjoint representation [jn|. For simplicity we 



concentrate on the SU(2) case. Then the generalization of (|9]) looks as: 

m,F, v {B) = m,*{n a F; v {A)) (10) 

where F^ u is the non-Abelian field strength tensor and n a is an arbitrary vector in the 
color space. Again, the choice of n a is a matter of gauge fixing and the physical results 
do not depend on this choice. One can visualize the vector n a as the direction in the 
color space of the magnetic fields transported along the Dirac strings attached to the 
monopoles. 

One can derive a Zwanziger-type Lagrangian which ensures the validity of the con- 
straint flltiD , see Section [4.3| . This formulation allows then to make predictions for the 
heavy monopole potential (|7|). In particular, at short distances the potential can be 
derived from the first principles: 

V mm {R) = -1-^, Qm-9 = 2n (11) 

where g is the standard coupling entering the non-Abelian Lagrangian. To appreciate the 
meaning of flTT| ) it is instructive to compare it with the standard prediction for the heavy 
quark potential: 

Note that the overall coefficient in front of the potential, that is 3/4, reflects the non- 
Abelian nature of the quarks which belong to a fundamental representation. There is no 
such factor in case of the monopole potential (]TT|) . Thus, monopoles in SU{2) gluody- 
namics are Abelian! 

On the other hand, the non-Abelian nature of the theory gets manifested in the running 
of the coupling. In the both cases ( JE2| ) and (|TTJ) the quantum corrections result in g 2 — > 
g 2 (R)- Although the final result is similar, the derivation is actually different. Namely, 
in case of the radiative corrections to flT2| ) one deals of course with a usual field theory. 
In case of (|TT|), performing the regularization procedure one should address anew the role 
of the Dirac string attached to the heavy monopoles and check that this procedure is 
respecting the Dirac constraints. This is one of most subtle and exciting points and we 
are going to discuss it in details in Sections |4.2| , |4 It is worth emphasizing that there 
appearing first direct measurements of the V mi?l on the lattice JZIJ so that the prediction 



PD can be checked through measurements [13 



The Abelian nature of the monopoles poses another serious problem which is the 
theoretical understanding of the so called Casimir scaling [^TJ. First, let us remind the 
reader what the Casimir scaling is. One can measure the heavy quark potential not only 
for the sources in the fundamental representation, as we discussed so far, but in any other 
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representation as well. At very large distances the potentials are expected to depend 
crucially on the representation. For example, if the quarks have isospin T = 1 then at 
large distances the confining string connecting the quarks should break into two mesons 
and the potential, respectively, is expected to flatten out. However at all the distances 



available for the lattice measurements so far an approximate equation for Vq q holds [^2 



V T (R) « -T(T + l)—^- + T(T+l)aR, (13) 

7T K 

where the string tension a does not depend on the isospin T. The Casimir scaling, at 
least at first sight, makes Abelian dominance very questionable. Indeed, in, say, T=l 
representation one quark is neutral with respect to any singled out U(l) subgroup of the 
SU(2) and this quark should have escaped any interaction, in bald violation of QIB]). 



We will argue |H| that in the approach where the choice of the U(l) subgroup is a 



matter of the gauge fixing, see above, the problem of the Casimir scaling can be settled. 
Within this approach the heavy quark potential is: 

OL 

V t {R) « -T(T + 1)— ^ + T(T + l)a T (m H ,m v )R, (14) 
where mn, my are the parameters of the dual model. The Casimir scaling corresponds 



then to a particular limit, ma I my ^> 1. In Section [4.6| we will discuss this problem in 
detail, along with other related phenomenological issues. 

In Section ||] we show that the results of the numerical calculations in lattice gauge 
theories confirm the dual superconductor model of the gluodynamic vacuum. Moreover 
the recent study of the confining string structure and monopole structure support the 
theoretical ideas given in this review. 

The review can be subdivided into two parts: Sections |2| |3| contain the introductory 
remarks and main definitions, the advanced topics are discussed in Sections |] - 15]. 



2 Monopoles as Classical Solutions 

In this section, which is pure pedagogical in nature we summarize the facts known about 
the magnetic monopoles viewed as solutions to the classical equations of motion. We 
consider both electrodynamics and gluodynamics. 

2.1 Dirac Veto, Dirac Quantization Condition 

As is mentioned in the Introduction, there is a number of constraints imposed on the 
theory to ensure that the Dirac string does not produce any physical effect. First, there 
is the Dirac veto which forbids any direct interaction with the string. The condition is 
not trivial at all, in fact. Indeed, let us consider perturbation theory and the scattering 
of a charged particle off the string. In perturbative approach one uses the basis of plane 
waves. Which implies that the particles can be found at any space-point and their wave 
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functions overlap with the string. Which is in violation of the Dirac veto. Thus, a simple 
minded Born approximation can be misleading. Of course, this consideration alone does 
not yet allow to judge whether this overlap is significant numerically. But a little bit 
more involved estimates do indicate that it is important and the perturbation theory 
can and does bring wrong results. We shall examine these issues in much more detail in 



Sections [4.2| , |4.4| devoted to the radiative corrections and here just state that a consistent 
treatment of the Dirac string is always non-perturbative in nature because of the Dirac 
veto. This simple observation allows to resolve many theoretical puzzles. 

The best known constraint on the Dirac string seems to be the Dirac quantization 
condition which ensures absence of the Aharonov-Bohm effect for the electrons scattered 
off the string: 

Adx = QeQm = 2nk, (15) 



Qe<f> 
J st 



' string 

where Q e is the electric charge of the electron and k is an integer number. Note that 
as a result of the quantization condition the interaction of the monopoles is described 
by a large coupling constant if the electrons interact weakly and vice versa. Also, the 
interaction of an electron with the monopole is always governed by a coupling of order 
unit and is never weak. This is another reason for the use of non-perturbative methods 
in the theory of monopoles, which is of course always a challenge. 



2.2 Energy of the Dirac String 



Unlike the quantization condition (|T5j) , the energy of the Dirac string is rarely discussed 
in the literature. It is worth to emphasize, therefore, that naively the energy of the string 
is infinite in the ultraviolet. 

Indeed, let us estimate the energy of the string: 



-string 



We accounted here for the fact that the magnetic flux is quantized (see above) 

2irk 



(flux) = J H • ds 



Qe 



and tended the cross section of the string denoted by (Area) to zero at the end of the 
calculation. Thus, we substituted (Area)^ 1 by Afj V . 

The radial part of the magnetic field is also associated with an infinite energy: 

e rad ~ f(H Coul ) 2 d 3 r ~ - ~ A uv . (17) 
J r 

Note that this ultraviolet divergence is linear, i.e. weaker than the divergence due to the 
string, see Eq. (fL6|). 
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The infinite magnetic field of the string may have more subtle manifestations as well. 
Consider interaction of two magnetic monopoles with magnetic charge ±Q m placed at 
distance R from each other. Then, by the analogy with the case of two electric charges, 
we would like to have the following expression for the interaction energy: 

f tt CouItjCouI j 3 Qm ^ /io\ 

e int = J H t ti 2 dr=_ ^T^- 

Note, however, that if we substitute the sum of the radial and string fields for H 12 , then 
we have an extra term in the interaction energy: 

e mt = J (Hf °^Hf r + Hf r H^ oul ) d 3 r = +2^|i. (19) 

In other words, the account of the string field would flip the sign of the interaction 
energy! This contribution, although looks absolutely finite, is of course a manifestation 
of the singular nature of the string mag netic field, |H 5tr | ~ (Flux) /(Area). Note that 
the integral in (|T9|) does not depend on the shape of the string. 

To maintain the unphysical nature of the Dirac string we should use a regularization 
scheme which would allow to get rid of these singularities. 



2.3 Lattice Regularization in Abelian Theories 

Since the monopoles are naively having divergent energy (or action) in the ultraviolet, the 
regularization is a crucial issue. Moreover, we would like to follow the lattice formulation 
since the monopoles are observed on the lattice. And for a good reason, as we shall 
immediately see. 



Consider first the U(l) case. As is emphasized in Ref. |ID[ , the lattice formulation 
implies that a Dirac string which produces no Aharonov-Bohm scattering costs no action 
as well. The reason is very simple. The lattice action is written originally in terms of the 
contour integrals like ([15]) rather than field strength F^: 

S = Re exp{zQ e / A^} , (20) 

p Jdp 

where the sum is taken over all the plaquettes p. Thus, the condition ([H]) means absence, 
in the lattice formulation, of both the Aharonov-Bohm effect and the quadratic divergence 
flip]). Moreover, it is straightforward to see that the interference term (|T9| ) also vanishes. 

Later, we will also discuss Dirac strings which correspond to negative plaquettes in 
the lattice formulation. Such strings are associated with the fundamental monopoles with 
Q m = 1 in units of the Dirac quantization condition discussed above. The energy of such 
Dirac strings is infinite in the continuum limit, in agreement with the naive estimate ([16|) . 
The interference term (|T9"D, however, disappears in the lattice formulation in this case as 
well. Moreover, this observation is crucial to understand the heavy monopole potential 
in the continuum limit. 
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The radial field, H^ "' may also cause problems with infinite energy, see (|17]). The 
lattice regularization is not much specific in that case, however. The role of ro is simply 
played by the lattice spacing a. Thus, the probability to find a monopole on the lattice 
is suppressed by the action as: 

e~ s ~ exp(—const-Q~ 2 L/a), (21) 

where L is the length of the monopole trajectory, and the Q~ 2 factor appears because of 
the Dirac quantization condition (|I5|) which relates the magnetic charge Q m to the inverse 
electric charge. 

Although the Eq. fl2~T|), at first sight, rules out monopoles as physically significant 
excitations, the fate of the monopoles in the U(l) case depends in fact on the value of the 
charge Q e . The point is that the entropy factor, or the number of various trajectories with 
the same length L grows also exponentially with the length of the monopole trajectory: 

(Entropy) ~ exp(+const' ■ L/a) , (22) 

where the const' is a pure geometric factor. As a result for Q e ~ 1 there is a phase 
transition corresponding to the condensation of the monopoles. This phase transition, 
which is well studied on the lattice, is the first and striking example of importance of the 
UV regularization in the non-perturbative sector. 



2.4 Classification of Monopoles in non-Abelian Theories 

From now on, we will discuss monopoles in unbroken non-Abelian gauge theories, having in 
mind primarily gluodynamics, i.e. quantum chromodynamics without dynamical quarks. 
Moreover, for the sake of simplicity we will consider only the SU{2) gauge group. 

A natural starting point to consider monopoles in non-Abelian theories is their clas- 
sification. There are actually a few approaches to the monopole classification and it is 
important to realize both similarities and differences between them. 

The dynamical, or U(l) classification. Within this approach f24jl , one looks for 
monopole-like solutions of the classical Yang-Mills equations. Where by the "monopole- 
like" solutions one understands potentials which fall off as 1/r at large r, see Eq. (^|). The 
basic finding is that there are no specific non-Abelian solutions and all the monopoles can 
be viewed as Abelian-like embedded into the SU(2) group. Moreover, using the gauge 
invariance one can always choose the corresponding U(l) group as, say, the rotation 
group around the third direction in the color space. According to this classification, the 
monopoles are characterized by their charge with respect to a U(l) which are integer 
numbers: 

\Q m \ = 0, 1, 2, ... . (23) 
The topological, or Z 2 classification. The Z 2 classification |25[ is based entirely on 



topological arguments. Namely, independent types of monopoles are enumerated by con- 
sidering the first homotopy group of the gauge group. The SU(2) gauge group is trivial 
since tci(SU(2)) = 0, while in case of 50(3), 

n 1 (SU(2)/Z 2 ) = Z 2 , (24) 



9 



and there exists a single non-trivial topological monopole. We will denote the magnetic 
charge of such monopoles as \Q m \ = 1. Note, however, that the charges Q m = ±1 are 
indistinguishable in fact. As for the charges Q m = 2 they are equivalent, from this point 
of view, to no magnetic charge at all. 

The topological classification (pi ) is readily understood as a classification of the Dirac 



strings whose end points represent the monopoles. Then there is only one non-trivial 
string, that is the one for which Eq. is satisfied for gluons but not for quarks. 

Namely, because the U(l) charge associated with gluons is twice as large as that of the 
quarks we may have 

exp{i# £ A^dx 1 *} = - 1 (25) 

and such a string is still not visible for the isospin-one particles. On the other hand, 
the standard plaquette action is based on the phase factor evaluated for particles in the 
fundamental representation. Which means, in turn, that the Dirac string is piercing the 
negative plaquettes. 

2.5 Z2 Monopoles 

In principle, the U(l) and Z 2 classifications are different. Indeed, while the U(l) classi- 
fication allows for any integer charge, the Z 2 classification leaves space only for a single 
non-trivial charge: 

Qm = 0,1. (26) 

The reconciliation of the two classifications is that the U(l) solutions with \Q m \ > 2 are 
in fact unstable because of the presence of massless charged vector particles (gluons) 0. 
The instability of the solutions implies that even if an external source with \Q m \ > 2 were 
introduced into the vacuum state of the gluodynamics, charged gluons would fall onto the 
monopole center because of the strong magnetic interactions. Moreover, one can imagine 
that as result of this instability the charged fields A are build up as well. 

In a somewhat related way, one can demonstrate the apparent irrelevance of the | Q m \ = 
2 monopoles by producing an explicit non-Abelian field configuration which looks as a 
\Q m \ — 2 monopole in its Abelian part but has no SU{2) action at all |13| . This field 
configuration is a Dirac string with open ends, which correspond to a monopole-anti- 
monopole pair separated by distance R. Note that the Abelian flux is still transported 
along the Dirac string and is still conserved for the radial field. What is lost, however, is 
the relation between the Abelian flux and action. In the Abelian case non-vanishing flux 
means non-vanishing magnetic field and non-vanishing action since the action density is 
simply H 2 . Now the action is {F^ u ) 2 and the Abelian part of the F? v can be canceled by 
the commutator term. 

The simplest example illustrating the cancellation between the Abelian component 
and the commutator term in F® v is produced by the potential: 

A a = _ , (27) 
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where a a are the Pauli matrices and r a is the radius vector. One can readily check by 
a direct calculation that the corresponding non-Abelian field-strength tensor vanishes 
identically (some care should be exercised to analyze the singularity at r = 0). Another 
way to convince oneself that the configuration (p7[) has zero action is to observe that 

a; = (28) 

where 

il° = ia a n\ 

and n a is the unit vector looking from the center of the "monopole" to the observation 
point. 

We pause here to note that the potential ([2j]) does not exhibit any Dirac string. The 
reason is that we have not gauge rotated it to the form suitable for the U(l) classification 
of the monopoles. The absence of the Dirac string may look appealing at first sight. On 
the other hand, a shortcoming of using an explicitly non-Abelian field, like (^) is that one 
cannot add up in this case the monopole fields if the monopoles are situated at different 
space points. Thus, we are not using such forms in this review although one cannot rule 
out that in future some progress can be made to overcome the difficulty with adding up 
the monopoles written in the form similar to ( f2~7| ) . 

Coming back to the problem of writing a filed configuration which in its Abelian part 
looks as a monopole-antimonopole pair (\Q m \ = 2) we will use formulation with a Dirac 
string. Algebraically the problem becomes more complicated and here we reproduce only 
the final answer Namely, such a configuration is generated from the vacuum by the 
following gauge rotation matrix: 



n=UVf^^J' (29) 

where <p is the angle of rotation around the axis connecting the monopoles and Ad is the 
£7(1) potential representing pure Abelian monopole pair: 

Apdxp = \ [^~~~ S j di P = A D (z,p)d(p, (30) 

where z± = z ± R/2, p 2 = x 2 + y 2 , r± = z±+ p 2 . Note that the action associated with 
the Dirac string is considered in this case zero, in accordance with the lattice version of 
the theory (for details see |T3| ). 



In this example, the monopoles with \Q m \ = 2 are a kind of a pure gauge field config- 
urations carrying no action. 

It is somewhat more difficult to visualize dynamically the equivalence of the Q m = ±1 
monopoles, also implied by the Z 2 classification. The mechanism mixing the Q m = ±1 
solutions seems to be the following. Imagine that we start with, say, Q m = +1 solution. 
Then a Dirac string carrying the flux corresponding to the Q m = —2 can be superimposed 
on this solution. It is important at this point that such a Dirac string costs no action (or 
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energy). Then the radial magnetic field can also change its direction since this does not 
contradict the flux conservation any longer. In a related language, one could say that the 
\Q m \ — 2 monopoles are condensed in the vacuum and as a result the magnetic charge 
can be changed freely by two units. 

As far as interaction of two \Q m \ = 1 monopoles is concerned, one might expect 
that they would interact as a monopole-antimonopole pair rather than as a monopole- 
monopole configuration. Indeed, the monopole and antimonopole attract each other and 
thus represent the lowest energy state of the system. 



2.6 Conclusions # 1 

Thus, the physics of the monopoles in non-Abelian theories in the classical approximation 
turns very simple. 

Namely, there exist only monopoles with \Q m \ — 1 = 27r/<? where g is the coupling 
constant of the non-Abelian SU (2) theory. The monopoles are infinitely heavy and their 
interaction is Abelian like: 

V - - - Vm - — cm 

""" ~ 1 rr/,' ~ g*R' ' 

where R is the separation between the monopoles. 

Clearly enough, this first approximation falls far beyond an adequate description of 
the empirical data on the monopoles, see the Introduction. Thus, we are invited to go 
into more advanced approaches which we would try to introduce step by step. 



3 Lattice Monopoles, Introductory Remarks 

3.1 Dirac Strings on the Lattice, Non-Abelian Case 

The Dirac strings remain essentially Abelian objects in case of non-Abelian theories as 
well. A novel feature is that the field transported along the string can be arbitrarily 
oriented in the color space. Which is a kind of a new gauge freedom. 

In more detail, the general one-plaquette action of 577(2) lattice gauge theory (LGT) 
can be represented as: 

SiatiU) = (3 $P (l - \TrU[dp]) , (32) 

where (3 = 4/g 2 , g is the bare coupling, dp is the boundary of an elementary plaquette 
p, the sum is taken over all p, U[dp] is the ordered product of link variables XJ\ along dp. 
In particular, if Sp(x) = x then (|32|) is the standard Wilson action. The exponent of the 
lattice field strength tensor F p is defined in terms of U[dp\: 

U[dp] = = cos^|F p |]+zr a n a sin[^|F p |], (33) 
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where F = F a ■ r a /2, \F\ = \/F a F a and we define n a p = F£/\F P \ for \F P \ ^ 0, while n a p is 
an arbitrary unit vector for \F p \ = 0. 

The lattice action ( j3"!|) depends only on cos [||i^,|]. Therefore the action of the SU{2) 
LGT possesses not only the usual gauge symmetry, but allows also for the gauge trans- 
formations which shift the field strength tensor by Ank, \F p \ — > \F p \ + Airk, k £ Z: 

e ip * = exp{i\F p \n p } = exp{i(\F p \ + 4vr)n p } = exp{i(F^ + Aim a p )T a /2}. (34) 

Thus, the symmetry inherent to the lattice formulation can be represented as: 

F p a -> F; + 47rn a p , F p x n p = 0, n\ = 1. (35) 

The symmetry (|35"D is absent in the conventional continuum limit, f(F® u ) 2 d A x. If n p ^ 
then one can say that a Dirac string is piercing the plaquette. Similar to the Abelian 
case, such a Dirac string costs no action. 

Note that in the continuum limit n p becomes a singular two-dimensional structure 
which is representing the Dirac string world sheet. 



3.2 The 't Hooft Loop on the Lattice 

So far we discussed invisible Dirac strings. As we know, the Dirac strings corresponding to 
the fundamental monopoles with Q m = 1 have an infinite action, or energy. In the lattice 
formulation they correspond to the plaquettes with the phase factor exp{m). This infinite 
energy precludes the fundamental monopoles from being dynamical objects present in the 
QCD vacuum. However, they can be introduced as external objects. The infinite energy 
of the corresponding Dirac string then does not matter, the same as with the infinite self- 
energy of heavy quarks introduced via the Wilson loop. The fundamental monopoles can 
be introduced on the lattice via the 't Hooft loop [0. In a simplified way, the monopoles 
are visualized as end-points of the corresponding Dirac strings which in turn are defined 
as piercing negative plaquettes. The trick to introduce the negative plaquettes on the 
lattice is to formally change the sign of the coupling g 2 on a manifold of plaquettes. Then 
these plaquettes become negative in the limit \g 2 \ — > 0. 

Proceeding to more detailed definitions, the 't Hooft loop is formulated [18[] in terms 
of the action 

S(f3,-(3)=(3Y, TrU p -(3j2 TrU p , (36) 

p^ai peM 

where M is a manifold which is dual to the surface spanned on the monopole world-line 
j. Introducing the corresponding partition function, Z(/3,—/3) and considering a planar 
rectangular T x R, T 3> R contour j one can define 

T/ m\ 1 1 ~g) (V7\ 

V mrh (R) = --In-^-. (37) 

By the analogy with expectation value of the Wilson loop the quantity V mitl {R) is referred 
to as monopole-antimonopole, or heavy monopole potential. The potential V m fn{R), in a 
way, is the same fundamental quantity as the heavy-quark potential Vqq and its under- 
standing within the fundamental QCD would be of great importance. 
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3.3 Lattice Monopoles in the Abelian Projection of SU(2) Glu- 
odynamics 

The SU{2) gauge fields XJ\ on the lattice are defined by SU{2) matrices attached to the 
links /. These lattice fields are related to the continuum SU(2) fields A: U X:fl = e m9A ^ x \ 
here a is the lattice spacing. Under the gauge transformation, the field Ui transforms as 
= £7+ U x ^ ft x +fi, the matrices of the gauge transformation are attached to the sites x 
of the lattice. 

The Maximal Abelian gauge is defined on the lattice by the following condition ||26|| : 

maxR[U?], RlUil^TrfaUfasUt], l = {x,fi}. (38) 
S2 i 

This gauge condition corresponds to an Abelian gauge, since R is invariant under the 
gauge transformations defined by the matrices: 

6 e -ia( x) ) , a e [0, 2tt) . (39) 

Let us parameterize the link matrix U in the standard way 

U t =( COS ^ i01 . shl ^ e Z), (40) 
' V - sm ipie~ lXl cos (fie~ l(>l J ' y J 

where 9, x £ [— tt? +?t) and if G [0, tt). In this parameterization, 

R[Ui] = £>s2 W . (41) 
i 

Thus, the maximization of R corresponds to the maximization of the diagonal elements 
of the link matrix (fEJ). 

Under the £7(1) gauge transformations, (|39]), the components of the gauge field ([40|) 
are transformed as 

Therefore, in the MaA the gauge, the field 9 is the U(l) gauge field, the field x is the 
Abelian charge 2 vector matter field, the field <p is the non-charged vector matter field. 

A configuration of Abelian gauge fields 9\ can contain monopoles. The position of 
the monopoles is defined by the lattice analogue of the Gauss theorem. Consider the 
elementary three-dimensional cube C (Figure |l](a)) on the lattice. 

The Abelian magnetic flux H through the surface of the cube C is given by the formula 



Z7r P&dC 

where 9p is the magnetic field defined as follows. Consider the plaquette angle 9p = 
9\ + 9 2 — 9 3 — 9 A = dQ, the #j's are attached to the links i which form the boundary of 
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t H 




Figure 1: The magnetic flux through the boundary of the cube C (a) and lattice field 
strength tensor (b). 

the plaquette P, Figure [j](b). The definition of Op is Op = Op + 27rfc, where the integer k 
is such that — 7r < Op < n. The restriction of Op to the interval (— 7r,7r] is natural since 
the Abelian action for the compact fields Oi is a periodic function of Op. Equation (fI5|) 
is the lattice analogue of the continuum formula m = § H dS. Due to the compactness 
of the lattice field (— n < < n) there exist singularities (Dirac strings), and therefore, 
divH ^ 0. 

The magnetic charge m defined by eq. (^) has the following properties: 

1. m is quantized: m = 0, ±1, ±2; 

2. If m ^ 0, then there exists a magnetic current j. This current is attached to the 
link dual to cube C. 

3. Monopole currents j are conserved: 5j = 0, the currents form closed loops on the 
AD lattice. 

We discuss in details the properties of the lattice monopoles in gluodynamics in Sec- 
tion 

3.4 Dirac Strings and Singular Gauge Transformations 

If we wish to build up the continuum theory in a way reproducing the basic features of the 
lattice formulation, we need first of all to ensure the invisibility of the Dirac string. One 
way is to put the energy of the Dirac strings to zero "by hands" . A more educated way 
is to allow for singular gauge transformations which would correspond to introduction of 
the Dirac strings. Let us start with the compact electrodynamics (cQED). 
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Under the gauge transformation, —>■ + c^a, the field strength tensor transforms 
as: F^y — > F^ u + [d^, d u ]a. For a singular gauge transformation [d^, d u ]a ^ 0. Moreover, 
for a sufficiently general singular gauge parameters a the commutator of derivatives is 
proportional to the 2D 5 function on some closed surface: [d^, d u ]a oc If x are 

the 4D coordinates of the surface, a a (a = 1,2) are 2D coordinates on the world-sheet, 
then the general representation of the surface is: 



J flU 



d 2 o^fg V(cr) <5 (4) (x - x{a)) , g(cr) = Det[ d a x^ dpx^ ] , 



(44) 
(45) 



Now, we can write down a partition function of electrodynamics which is invariant under 
the singular gauge transformations: 



(46) 



Of course / PS is only formal notation for the summation over all surfaces in 4D as far 
as we do not specify the measure. On the other hand, this summation is well defined on 
the lattice and fljE|) can be viewed as the partition function of cQED with the Villain type 
of the action written in the continuous notations, if X>£ in fl4~6|) includes summation over 
closed and not closed surfaces. 

Proceeding now to the non-Abelian case, we note again that the symmetry (F~ 



V is 

absent in the conventional continuum action, J(F^ U ) 2 d A x. To bring the continuum the- 
ory into agreement with the lattice formulation, we need a generalization of fl46|) . In 
the continuum limit becomes a singular two-dimensional structure = -^e^xpT^ 
which is a generalization of the Dirac strings in the compact electrodynamics and which 
transforms in the adjoint representation of the gauge group: 



d 2 a^-gt a (a) 6®(x-x(a)) 



(47) 



The surface (f47D need not to be closed. Moreover, the second equality in ( |35|) requires 
that 

V,(a) x *F^(x) =0, (48) 

where the continuum field strength tensor F^ v = d^A v — d v A^ — i[A^A v ]. Eq. (f4"5|) 
determines the color structure of : 



V((j) n a {a) 



n a {a) = (t-*F a ) (t-*F 



-1/2 



(49) 



where (t ■ F a ) = t^ v {a) F^ u {x) and n a is normalized as n 2 = 1. On the set of points where 
(t ■ *F a ) = the direction of n a (a) is arbitrary. 

Now, the continuum analog of the lattice symmetry (|35|) is: 



F a + 47r *S a 



(50) 
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The action of SU(2) gluodynamics which possesses the additional symmetry fl50f) can be 
formally represented as: 

Z = JvA exp{-£(F)} , (51) 
' r ^ (52) 



S(F) = - log J PSexp{-— J d A x 
where the integration is over all possible surfaces, 



'} 



d 2 a iiv n a (a) 5 {4 \x - x{a)) 



(53) 



The expressions (^Tj|52 



are only formal since it is impossible to separate rigorously 
the measure from the gauge degrees of freedom in DA. Nevertheless, the Eq. fl5T| , |5"2] ) 
is a good starting point since it reproduces the basic symmetries of the theory. Indeed, 
the action is invariant under smooth £77(2) gauge transformations since vector n a 
transforms in the same way as does. It is also invariant under transformations (|50|) 
which correspond to the lattice symmetry relations (|35|). 

Note also that for self- intersecting surface Eq. (]55|) , the world-sheet vector field 
n a (a) is generally multi- valued as function of x. Furthermore, for the non-orientable 
surfaces the field n a (a) cannot be defined smoothly everywhere on S. To avoid these 
complications we consider only the orientable surfaces without self-intersections. This 
reservation is specific for Dirac strings in the non-Abelian case. 

3.5 The 't Hooft Loop in the Continuum Limit 

After the experience with working out the continuum limit of the Dirac string world sheet 
in the preceding subsection, there is no difficulty to figure out that a continuum analog 
of the 't Hooft loop looks as [ |13]] : 



exp 



40= 



d x 



pa 



^+27^; 2 



(54) 



d n° 



6 {4 \x-x{a)) 



(55) 



where the surface spanned on the contour C is assumed to be non-intersecting. The 
unit three-dimensional vector field n a (a), n 2 = 1 is defined on the world-sheet of the 
Dirac string. For analogous definitions see (f49|). Therefore n a (a) is not an independent 
variable, it is completely determined by the components of the field strength tensor F£ v . 
It can be shown |13| that the Eq. (p4|)-(^9|) define the correct 't Hooft loop operator the 
expectation value of which depends only on the contour C, not on a particular position of 
the surface S c . 
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3.6 Conclusions #2 



The main lesson from the lattice formulation which we learnt so far is that to match 
the lattice and continuum formulations of the gluodynamics one should allow for singular 
gauge transformations. In this way one can hope to reproduce the compactness of the 
U(l) subgroups of the SU(2) group and the non-observability of the Dirac strings which 
is crucial for introduction of the monopoles. 



4 Lagrangian Approach to the Dual Gluodynamics 

In this section we are going to derive some important results concerning monopoles and 
their interaction in gluodynamics, understood as a continuum limit of the lattice for- 
mulation. At first we discuss the Zwanziger receipt of introduction of the dual photons 



(Subsection |4.1|) . In Subsection [O^ we introduce a Zwanziger-type Lagrangian for gluody- 
namics. The central point here is the match of ordinary gluons in an adjoint representation 
with a single dual gluon which is an U(l) Abelian gauge boson. In Subsections |47|, |4.4| 
we use the Lagrangian formulation to derive the short distance behavior of the heavy 
monopole potential, with one-loop quantum corrections included. In Subsections [L5], fi6] 
we discuss the phenomenological Lagrangians describing the monopole condensation. The 
central point here is a derivation of the Casimir scaling. 

4.1 Photons, Direct and Dual 

Now, that we know that magnetic monopoles are relevant degrees of freedom on the lattice, 
at least in some gauges, we would like to have a continuum version of gluodynamics which 
would allow both for color and magnetic charges. We begin, however, with a review of 
the corresponding construction in case of the electrodynamics, that is of the Zwanziger 
Lagrangian ||T^| . 



We have already mentioned some basic features of the Zwanziger approach in the 
Introduction and here would like to directly proceed to the Lagrangian: 

L Zw (A,B) = ^(m-[dAA}) 2 + i(m • [d A B]f + (56) 

+ Um- [d A A])(m ■ *[d A B}) - Um ■ [d A B])(m ■ *[d A A]) + ij e -A + ij m -B, 

where j e ,j m are electric and magnetic currents, respectively, m p is a constant vector, 
m 2 = 1 and 



[A A B\ LV = A p B v - A V B^ , (m-[AA% = m v [A A B] 

*[A A B\^ v = ^s, uXp [AAB] Xp . 

At first sight, we have introduced two different vector fields, A and B, to describe inter- 
action with electric and magnetic charges, respectively. If it were so, however, we would 
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have solved a wrong problem because we need to have a single photon interacting both 
with electric and magnetic charges. And this is what is achieved by the construct (|56|). 
Indeed, the action ( |56l) is not diagonal in the A, B fields and one can convince oneself 
that the form of the A, B interference terms in fl5"6|) is such that the field strength tensors 
constructed on the potentials A and B are in fact related to each other: 

m^ u {A) = m*F^,{B) , (57) 



where *F denotes the dual tensor. In perturbation theory, Eq. (|57j) ensures that there are 
only two physical degrees of freedom corresponding to the transverse photons which can 
be described either in terms of the potential A or B. Topological excitations, however, 
can be different in terms of A and B. 

The physical content of ([56]) is revealed by the propagators for the fields A, B. In the 
a-gauge one can derive: 

(A,A V ) = (B,B U ) = 1 (<V - (1 - a) , (58) 

(59) 

i 

k 2 (km) 



(A,B U ) = -{B,A U ) = ^>A^. (60) 



The propagators should reproduce, as usual, the classical solutions. And indeed, the 
(AA), (BB) propagators describe the Coulomb-like interaction of two electric charges and 
magnetic monopoles, respectively. While the (AB) propagator reproduces interaction of 
the magnetic field of a monopole with a moving electric charge. The appearance of the 
poles in (k ■ m) is a manifestation of the Dirac strings. 



To summarize, the Zwanziger Lagrangian in electrodynamics |T9[ reproduces the clas- 
sical interaction of monopoles and charges. Upon the quantization, it describes the correct 
number of the degrees of freedom associated with the photon. 

4.2 Monopole Charge Renormalization in Abelian Theory 

We will now demonstrate that the quantum field theory with the Zwanziger Lagrangian is 
not self consistent, and the knowledge of propagators (^) is not enough to get the correct 
answers. To demonstrate these difficulties we consider the simplest theory containing 
light electrons, photons, and magnetically charged particles which are heavy and thus we 
can neglect their contributions to the vacuum polarization. 

Consider the renormalization of the monopole charge due to the vacuum polarization 
induced by the electron loop, see the review |§ and references therein. As is emphasized in 
Ref. |13|], the crucial point is that only loops with insertion of two external (i.e., monopole) 



fields can be considered despite of the fact that there is no perturbative expansion in g ■ e. 
Indeed, considering more insertions makes the graphs infrared sensitive, with no possibility 
for In A uv to emerge. 

Then, the evaluation of, say, first radiative correction to the propagator (B^B V ) in the 
Zwanziger formalism (p8|) seems very straightforward and reduces to taking a product of 
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Figure 2: (a) Renormalization of the gauge field propagator, (b) Renormalization of 
the gauge field propagator. 



two (AB) propagators and inserting in between the standard polarization operator of two 
electromagnetic currents, see Fig. 1(a). The result is ||: 

(B,B v ){k) = 6 -^(l-L) + —l—Qr-mpmjL, (61) 

where 

L = ^\nA 2 uv /k 2 

and we have neglected the electron mass so that the infrared cut-off is provided, in the 
logarithmic approximation, by the momentum k. 

At first sight, there is nothing disturbing about the result fl6T|). Indeed, we have 
a renormalization of the original propagator which is to be absorbed into the running 
coupling, and a new structure with the factor (k ■ m)~ 2 in front which is non- vanishing, 
however, only on the Dirac string. The latter term would correspond to the renormaliza- 
tion of the Dirac-string self-energy which we do not follow in any case since it is included 
into self-energy of the external monopoles. What is, actually, disturbing is that according 
to QBTD the magnetic coupling would run exactly the same as the electric charge (see Fig. 
1(b)), 

(A,A v )(k) = (62) 

violating the Dirac quantization condition. 

The origin of the trouble is not difficult to figure out. Indeed, using the propagator 
(AB) while evaluating the radiative corrections is equivalent, of course, to using the full 
potential corresponding to the Dirac monopole Ap. Then, switching on the interaction 
with electrons would bring terms like Af, ip'jip . Since A^ includes the potential of the 
string, electrons do interact with the Dirac sheet and we are violating the Dirac veto 
which forbids any direct interaction with the string. 

Let us demonstrate quantitatively that, indeed, it is the interaction with the Dirac 
string that changes the sign of the radiative correction. This can be done in fact in an 
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amusingly simple way. First, let us note that it is much simpler to remove the string 
if one works in terms of the field strength tensor, not the potential. Indeed, we have 
H = Ji Coul 4- H str while in terms of the potential any separation of the string would 
be ambiguous. 

We discussed already that the account of the string field would flip the sign of the 
interaction energy, see eq. (|19|) . On the other hand the duality principle ([OJ) implies that 
at the classical level we should write: 

Hx-H 2 = Hf ouZ -H^. (63) 

However, the electrons in the diagram shown in Fig. 1(a) interact with the full po- 
tential and therefore the first radiative correction would bring the product of the 
total Hi • H2 which includes also the string contribution^. Indeed, the result in the log 
approximation would be as follows: 

<J(Hi • H 2 ) = L(Hf r + H.c° ul ) ■ (Hf r + H.% oul ) = - LHF x oul ■ , (64) 

where at the last step we have used the observation (|T9|). 

Now, it is clear how we could ameliorate the situation. Namely, to keep the Dirac 
string unphysical we should remove the string field from the expression (|64]) which arises 
automatically if we use the propagators (pq) following from the Zwanziger Lagrangian. 
Thus, one should change Hi • H 2 in the expression fl6"4|) to Hf oul ■ Hf? "'. The justification 
is that we should remove the effect of the string field from any observable and consider 
this as a constraint on the ultraviolet regularization. The constraint is satisfied by the 
lattice regularization. The corrected propagator [B^B V ) at the one-loop level is becoming 
then: 

(BpBvKk) = ^(1+L). (65) 

In particular, there is no Dirac string self-energy renormalization. 

The most important conclusion is that fl65|) does satisfy the Dirac quantization con- 
dition also for the running coupling (compare eq. ( |62|) and eq. (|65|) ). 

It might worth emphasizing that the running of the magnetic coupling became possible 
only because we do not have any longer the Maxwell equation div H = on the operator 
level. Indeed, we have a source of the dual photon (or gluon) which is mixed up with the 
ordinary photon (gluon). Similarly, for the electric field we have div E = p e \. Classically, 
for a point-like charge p e i = Q e S^(r) and the coupling is not running of course. However, 
on the quantum level the fact that div~E 7^ as far as the operator equations are concerned 
gets revealed through the running of the electrical coupling. Now the story repeats itself, 
in the dual formulation, for the magnetic coupling. In Subsection |4.4| we demonstrate 
that in the non-Abelian case the Dirac quantization condition is also valid for the running 
coupling. 

1 At this point we assume in fact that Ajjv is larger than the inverse size of the string, which is 
convenient for our purposes. Other limiting procedures could be considered as well, however. 
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4.3 Dual Gluon as an Abelian Gauge Boson 

Now, if we try to extend the Lagrangian approach to the dual gluons, we immediately 
come to a paradoxical conclusion that the dual field, if any, is Abelian. Indeed, we have 
already emphasized that monopoles associated with, say, SU(N) gauge group are classified 
according to [/(l)^" 1 subgroups |24| and can be realized as pure Abelian objects. Thus, 



there is no place for a non-Abelian dual gluon because the monopoles do not constitute 
representations of the non-Abelian group. 

From this simple observation we immediately conclude that the dual field B is an 
Abelian gauge field. Therefore in the Lagrangian a la Zwanziger the dual field B can be 
mixed up only with one of the fields A a . Moreover, which field A a is being mixed up is 
matter of a gauge fixing. The Lagrangian realizing these ideas looks as (for more detail 
see 11): 

1,™. , 2 1/ r „ „, „X2 



L dual (A a ,B) = -{F; v y + -{m-[dAB} - i*F a n a ) + % Jm B + i f e A a . (66) 



Now, a few remarks concerning the Lagrangian ( pq ) are in order: 

(a) First, if the magnetic current is vanishing, j m = then the integration over the 
field B reproduces the standard Lagrangian of the gluodynamics (we are omitting the 
gauge- fixing and ghost terms). 

(b) As far as the quantization is concerned, the Lagrangian (|66D reproduces the correct 
degrees of freedom of the free gluons. Indeed, in the limit g — > and for n a = 5 a,s the 
Lagrangian (|66| ) becomes: 

L dual (A a ,B) = l(d A A 1 ) 2 + ±{d A A 2 ) 2 + ij m B + if e A a + 

+ ^[m-(dAA 3 )} 2 + ^[m- (d A B)} 2 + (67) 
+ ^[m • (8 A A 3 )][m ■ *(d A B)} - % -[m ■ *{d A A 3 )}[m ■ {d A B)} , 



which is essentially the Zwanziger Lagrangian ( pq ) plus non-interacting non-diagonal glu- 
ons. Quantization at this point is the same as in the case of a single photon. 

(c) The emergence of the vector n a in the Lagrangian ( pq 1 ) is of crucial importance. 
The point is that the origin of the vector n a goes back to choosing the color orientation 
of the monopoles. As is emphasized above the monopole solutions are Abelian in nature 
which means, in particular, that they can be rotated to any direction in the color space 
by gauge transformations. Thus, picking up a particular n a is nothing else but using the 
gauge fixing freedom. Therefore, we can either average over the directions of n a or fix n a 
but evaluate only gauge invariant quantities, like the Wilson loop. 

(d) The Z 2 nature of the monopoles is manifested in the freedom of changing n a — > 
— n a , B^ — > —Bp. Indeed, under such transformation a monopole with the charge Q m = 
+1 is transformed into a monopole with Q m = — 1 and vice versa. In the language we 
used above such a transformation corresponds to adding a Dirac string with a double 
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magnetic flux. We see that the averaging over ±n a is a part of the overall averaging over 
all possible embedding of the U(l) into the SU(2) gauge group. 

(e) In its generality, the Lagrangian (|66|) obviously possesses SU(2) x £7(1) gauge 
invariance: 

517(2) : " f^f , ^(1) = B, -> B, + d,a . (68) 

All the features (a)-(e) indicate that the Lagrangian ( |66|) is a reasonable choice to 
describe the dual gluodynamics. Now, we would like to put the construction (|66"D on 
further tests. An obvious check is the evaluation of the quantum corrections to the 



heavy monopole potential V^ m at short distances. Indeed, it can be shown, [13], that the 
radiative corrections result in the standard running of the non-Abelian constant g 2 . 



4.4 The 't Hooft Loop and the Dual Gluon 

We have already derived a continuum analog of the 't Hooft loop, see Eq fl54|) . Now, using 
the notion of the dual gluon field we will turn to a more detailed analysis of the heavy 
monopole potential. 

For a rectangular contour T x R we can write for the surface S a : 

E« = n a -^—[mAj] (69) 
J (mo) 

which is a particular solution to <9£ = j. It is obvious that the current j will serve now 
as a source of the dual gluon B. Upon substitution of (^) the expectation value H(j) is 
represented as: 

H(j) = \\VA Jvn exp{- i ^|d 4 x(F a + 27rn a ^-*[mAj]) 2 }, (70) 



n 2 =l 



Note that Eq. (^) extends the definition of the field n a into the entire space-time 
while only its value of n a on the string Sj does matter, while a particular way of extending 



n a (a) — > n a (x) is irrelevant. Moreover, the path integral in (|70|) is to be performed with 
a constraint 

e abc n b (F c ■ j-^\*[m A j]) = 0, a = 1,2, 3, (71) 
(mo) 

which is a consequence of the constraint (49). The constraint equation fl7lf) may be 
implemented by an additional field x a '- 

H(j) = | J VnVAVx exp{ ~ J d 4 x (F a + (72) 



n 2 =l 



-^n a -^—*{mAj}) 2 + t—[jC} 
(mo) g J 
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C, = C,( X ,n,F) = —e^ x * n \m-*F c ),. (73) 

It is worth emphasizing that the field n a (x) is a kind of a fake variable in the representation 
(|T2"D. The path integral is clearly independent on n a away from the string fl69|), but at 
the same time n a (x) for x e £j is determined through ([71]). 

Coming back to the evaluation of the potential V^ m , on the classical level we tend the 
coupling g — > 0. In this limit and for n a = 5 a ' 3 the Lagrangian fl66|) reduces to: 

L° = hdAA°f + Um-[dAB - i*d A A 3 ]) 2 , (74) 

which essentially coincides with that of Zwanziger (0). In particular, to the lowest order 
in g 2 the monopole-antimonopole interaction is given by a single 5-field exchange. Using 



Eq. (58) one immediately obtains: 

1 ,2ti\o r . . i . , , m 4vr 2 



H(j) = exp{ -- (-) 2 / j,A-% } ~ exp{ T—A^)(R) } (75) 
^ 9 J 9 



Vmm(R) = ~~ 2 -^ (76) 
9 R 



Which is the expected result of course. 

Having full Lagrangian (|66D we can consider the radiative corrections to (|76|). The 
calculations similar to that of Section |4.2j give the result [13|] : 



(77) 

which is valid of course at small distances. Eq. (|77|) makes manifest that monopoles in 
gluodynamics unify Abelian and non-Abelian features. Namely, the overall coefficient, 
ir/g 2 is the same as in Abelian theory, while the running of the coupling g 2 {R) is the same 
as in the non-Abelian theory. The U(l) normalization is based on the classification of the 
monopoles in non-Abelian theory (see, e.g., and references therein). The running of 
the coupling reflects the general rule that the effect of the fluctuations at short distances 
can be absorbed into the renormalization of the coupling. Eq. (JTf\) also show that the 
Dirac quantization condition is valid at the perturbative level in non-Abelian theory. 

In order to derive flT7| ) we have to neglect the contribution of the Dirac strings. We 
should subtract the contribution of the strings, even if it is not vanishing within pertur- 
bation theory. Essentially, the rule is that the product ■ H 2 should not include the 
piece H. Coul ■ B_ str . 

To justify this rule we have to go actually beyond the Lagrangian (|66|) and consider 
a regularization procedure which takes into account the Dirac strings. It is worth em- 
phasizing that within the lattice regularization the term Ji Coul ■ H Sir is indeed absent. 
Let us recall the reader that in the lattice regularization the Dirac string pierces negative 
plaquettes. This is true in the limit g 2 — > 0, or a — > 0. If one looks for small deviation 
from U p = —1, then the action does not contain terms, which are linear in perturbations. 
Note that this will be not true for the expansion around arbitrary U p ^ ±1. This result 
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means in turn that in the continuum limit there is no term ~m Coul ■ H str . Thus, our naive 
removal of the effects of the Dirac string from the radiative corrections, is justified by the 
lattice regularization. 

The prediction of the monopole potential at large distances is also possible [[HJ if we 
use the Abelian dominance model. The existing numerical data |2(J do not contradict 
the predictions of Ref. 



T3|j. Higher statistics is needed, however, to crucially check the 



theory. 



4.5 Phenomenological Lagrangians 

We now switch to the physics of large distances, or the dual superconductor model of the 
confinement, see the Introduction. 

Here, the basic lesson brought about by the lattice simulations is that the monopoles 
condense (see reviews || and Section ^). In Subsection ^73] we have sketched theoretical 
arguments IJJ for the monopole condensation within the compact U(l). In case of the 



non-Abelian theory, a simple-minded generalization of these arguments is as follows. Let 
us begin with a lattice formulation and a small coupling g 2 SU (2)- Let us furthermore 
integrate, a la Wilson, over small-scale fluctuations and go in this way to a larger-size 
lattice. The only effect of this is the rescaling of the coupling g 2 su ^ 2 ) according to the 
rules of the renorm-group. However, we could expect that this procedure would not work 
further once we reach gf^) ~ 1- Indeed, the same non-Abelian coupling g"su(2) governs 
the physics associated with any U{\) subgroup. And we know that for gfj^i) ~ 1 there 
is a phase transition due to the monopole condensation. Thus, we expect the monopole 
condensation in non-Abelian theory as well, unless something else happens at smaller 

9su(2)- 

In the continuum limit and in the field theoretical language, one describes usually the 
monopole condensation within a (dual) Abelian Higgs model (AHM), see, e.g. reviews ]7[ 
|5|, U and references therein. The action of the model looks as: 

S= id i x{^Fl + \\{d-iB)^ + -\m 2 -V 2 f) (78) 



where Q m is the magnetic charge, is an Abelian gauge field, is the corresponding 
field-strength tensor, F^ u = d^B u — d u B^, and A, r\ are constants. The magnetically 
charged scalar field $ condenses in the vacuum, < $ >= rj, and the physical vector and 
scalar particles are massive, my = Q 2 r jf,mH — 2\rj 2 . 

The corresponding equations of motion possess an Abrikosov-type string solution. An 
open string with electrically charged particles at the end points is thought to mimic a 
confined QQ pair in QCD. The string tension is a function of the parameters of the theory, 
that is my, ttih- 

There is one subtle point about the dual superconductor model which we would like to 
emphasize. Namely, if we turn back to the simple arguments of the Subsection [2.3| which 
explain the monopole condensation, we can easily visualize that, indeed, the monopole 
condensate is non- vanishing, that is<$>=?77^0. Moreover, it is then consequence of 
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the gauge invariant interaction that m v ^ as well. However, it is not easy to understand 
at all which parameter is related to the Higgs mass, ma- One allows usually for arbitrary 
my,mH which are treated as fit parameters to reproduce the lattice data. Moreover, 
the most recent fits to the transverse structure of the confining string are close to the 
Bogomolny limit, my = m H see, e.g., [^jj and Subsection |5]2|. Of course, there is an 
intrinsic uncertainty in all the fitting procedure since one using the classical solutions 
to ( f78|) to fit the results of the lattice simulations which account fully for the quantum 
effects. 

From purely theoretical point of view another fitting procedure would be more logical. 
Namely, at large distances it is logical to assume the London limit since the monopole 
condensation in the compact U(l) does correspond to this limit. As for the finite Higgs 
mass itlh it should be produced by account for the non-Abelian gluons which are to 
be manifested at short distances. However, there is no theoretical papers which would 
produce an explicit procedure of this type and we cannot include any detailed discussion 
of such an approach into this review. Although, as we shall see in the next Subsection, 
further exploration of this approach seems urgent. 



4.6 Casimir Scaling 

As we already mentioned in the Introduction, a direct use of the Abelian models (|78| ) is in 
contradiction with the Casimir scaling. Indeed, whatever U{1) subgroup we choose there 
are quarks neutral with respect to this subgroup once the (color) isospin of the quarks is 
integer. However, after the experience with evaluation of the heavy monopole potential 
(see above) this puzzle does not look very difficult to solve. Namely, we can use a U(l) 
language for dual gluons. However, we should understand the choice of a particular U(l) 
as a gauge fixing procedure and evaluate gauge-invariant quantities. 

Thus, combining the idea of introducing an effective Higgs field $ with fundamental 
Lagrangian (|66|) we come to the following effective theory of SU(2) gluodynamics: 

Z = J VAVBV<5> V X J Vn exp{ - J d 4 xL eff } , (79) 

n 2 =l 



^cff 



-(F a 
4 V ^ 



+ 2 



m ■ [d A B 



>pa n a _ £ abc x a n b* F c 



+ 



SO) 



1 Air 

+ \(d+t-B)^\ 2 + vm), 

2 g 



where 



V(|$|) = A(|$| 



n 



and A and r\ are phenomenological constants. Of course, the vacuum expectation value 
of the Higgs, or monopole field is of order Aq C d- The physical assumption behind (fT9|) 



is that the effective size of the monopoles with \Q n 
although generically it is of order Aq C d (see also 



2 is in fact numerically small, 
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So far we discussed only embedding of the U(l) into SU(2). As far as the physics of 
large distances is concerned, there is a subtle point, how to define the U(l) dynamically. 
As is mentioned above the U(l) arises in lattice simulations through gauge fixing. The 
monopole properties depend on the particular choice of U(l). The Maximal Abelian 
projection turns most successful as far as the interaction of the heavy quarks in the 
fundamental representation is concerned ||. Thus, we will assume the use of this gauge. 
Moreover, let us fix the orientation of the dual U(l) group in the color space as the 
rotations around the third axis. Then knowing the classical, Abrikosov-type solution 
allows to find directly the string tension, <jAb{2T 3 ) as a function of the third component 
of the isospin of the quarks, T 3 and of the parameters of the Lagrangian([79|) . The role 
of the Wilson loop is to ensure that the QQ system is in the state with the total isospin 
T = 0. As a result, the string tension measured via the Wilson loop is given by: 

aT = 2TTI JV^ 27 ^' (82) 

where a Ah denotes the Abrikosov string tension, calculated in pure Abelian Higgs model 
with and for external charges of the magnitude ±2T 3 while the overall factor (2T + is 
due to the normalization of the wave function. Note that quantum mechanically Eq (|52]) 
is similar, say, to expressing the hyperfine splitting for a state with a given total spin in 
terms of the spin-spin interaction of the constituents. 

In the London limit, m# 3> my, the string tension UAb{2T 3 ) ~ Tf and we reproduce 
the Casimir scaling: 

1 T 1 

a T = (const) 2T - J2 T 3 = (const) - T(T + 1) . 

T 3 =—T 

In the Bogomolny limit, rrin = my, aAb(2T s ) ~ |T 3 | and o~t ~ T(T + 1)/(2T + 1). 

As is mentioned above, the structure of the confining string for the quarks in the fun- 
damental representation is best described in the Bogomolny limit, see Subsection |5]^ and 
Ref. [p7| . As is mentioned in the preceding Subsection, a possible way out of the difficulty 



could be account for the non- Abelian gluons at short distances. Indeed, the string tension 
is sensitive to much larger distances than the string structure in the transverse direction. 
At large distances, where the picture of the monopole condensation applies, one expects 
the London limit to hold (see the preceding Subsection). At smaller distances the charged 
gluons could become important. 

We conclude this subsection with a remark, that from the theoretical point of view, 
the most important limitation in the use of ( jSOf) is that it can be consistently treated on 
the classical level only. The difficulty to extend it to the quantum level is due to the Dirac 
strings. Indeed, perturbatively the Dirac veto is violated for virtual particles (see Subsec- 
tions |4.2j , O). When the monopoles condense the Dirac strings are filling the whole of the 



vacuum and there are no known ways to rectify the perturbation theory. In case of the 
Abelian Higgs model, it is even more convenient to use the dual language when the Dirac 
strings are attached to the electric charges. Then, to respect the Dirac veto one should 
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impose the condition that the Higgs field vanishes along a line connecting charges . 
The static Abrikosov-Nielsen-Olesen string satisfies this constraint ||. However, there is 
no known complete set of solutions satisfying this boundary condition. 



4.7 Conclusions # 3 

The use of the Zwanziger Lagrangian is a convenient means to describe interaction of 
photon with both magnetic and electric charges. Formally, one introduces two potentials, 
A^, Bp. However, because of the condition F^A) = *F^ U (B) there is a single photon in 
fact. One can get rid of one of the potentials. But then the Dirac string emerges explicitly 
and the formulation becomes non-local. 

Another point, very important from the point of view of the phenomeno logical ap- 
plications, is that in case of the gluodynamics the dual gluon (that is the field B) is 
still Abelian like although the "direct" gluons (that is the fields A a ) are in an adjoint 
representation of the non- Abelian group. 

An important conclusion concerning the phenomenological Lagrangians is that obser- 
vation of the Casimir scaling imposes the London limit, as far as one fits the lattice data 
with classical solutions. To avoid contradictions with fitting the structure of the confin- 
ing string one should assume then that the applicability of the AHM is limited at short 
distances by existence of the non-Abelian gluons, while the parameterization in terms of 
a finite Higgs mass, m# can be misleading. However, at present time there exists no 
explicit realization of this idea. 



5 Monopoles in Lattice Gluodynamics 

In this section we give the brief review of the numerical results obtained for monopoles 
on the lattice. We pay the special attention to two new topics: structure of the confining 
string and structure of monopoles in the Abelian projection. Most of the results are 
obtained for 577(2) lattice gauge theory for the maximal Abelian projection. The main 
definitions are given in Section |3.3| . 



5.1 Properties of Abelian Monopoles in SU(2) Lattice Gluody- 
namics 

At first we list the main properties of the Abelian monopoles in SU(2) lattice gluodynam- 
ics in the maximal Abelian projection ||26|| , the reader can find many important details in 



reviews [|9||. 

1. In the confinement phase the monopole currents form a dense cluster, we call it in- 
frared (IR) cluster, and there is a number of small mutually disjoint clusters, "ultraviolet" 
(UV) clusters. In the deconfinement phase the monopole currents are dilute. In Figure |3] 
we demonstrate these facts showing the abelian monopole currents for the confinement 
(a) and the deconfinement (b) phases. The IR cluster percolates [^0| and has a nontrivial 
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fractal dimension 
of the IR cluster, 
confinement of quarks |32 



3lj , Df > 1. The properties of UV clusters differs much from those 



it can be shown that the IR monopole cluster is responsible for the 





Figure 3: The abelian monopole currents for the confinement (a) {(3 = 2.4, 10 4 lattice) 
and the deconfinement (b) phases {(3 = 2.8, 12 3 • 4 lattice). All monopole currents are 
closed (conserved) due to the periodic boundary conditions. 



2. The SU(2) string tension is well reproduced by the contribution of the abelian fields 
and/or abelian monopole currents ||12j| . This property is called monopole (Abelian) dom- 
inance. While the Abelian dominance is an almost trivial property of the system ||33|| , the 
monopole dominance indicates that Abelian monopoles are important degrees of freedom 
for the infrared physics. 

3. The monopole condensate is nonzero in the confinement phase of SU(2) gluody- 
namics and it vanishes at the critical temperature corresponding to the deconfinement 
phase transition, thus it plays the role of the order parameter [[34j| . We can treat this fact 
as a justification of the dual superconductor model of the QCD vacuum, the monopoles 
playing the (dual) role of the Cooper pairs. 

4- The monopoles are correlated with the density of SU(2) action. The total action 



of SU(2) fields is correlated with the total length of the monopole currents [35], so there 
exists a global correlation. The Abelian monopoles in the MaA projection are also locally 
correlated with the non- Abelian action density |3^| . This fact shows that monopoles are 
some physical objects (not the artifacts of the singular gauge transformation), since by 
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0.03 




Figure 4: The fit of the lattice data by the classical solution of the discretized version of 
the Abelian Higgs model (ff8|). 



definition we call the object physical if it carries the action. 

5. The correlations of the monopole currents, the electric currents, the topological 



charge density and the action density was found in Refs. [37 



6. In Refs. [j38f the effective Lagrangian for monopoles was reconstructed from numeri- 



cal data for monopole currents for SU(2) gluodynamics in the Maximal Abelian gauge. It 
occurs that this Lagrangian corresponds to the Abelian Higgs model, the monopoles are 
condensed and the classical string tension of the Abrikosov-Nielsen-Olesen string describes 
well the quantum string tension of SU (2) gluodynamics. It means that the monopole de- 
grees of freedom are important for the description of the gluodynamics at large distances. 

5.2 Anatomy of the Confining String 



The authors of Ref. |39| presented the results of the numerical study of the confining string 
in SU (2) lattice gluodynamics. They measured the expectation value of the electric field, 

— t —* 

E, and the expectation value of the monopole current, k, near the line connecting the test 
quark-antiquark pair. For the test quarks placed at the z axis, only the z component, 
E z (p), of E and the azimuthal component, ^e(p), of k were found to be non- vanishing 
(here p = \/x 2 + y 2 is the distance from the center of the flux tube and the azimuthal 
angle 6 is defined as usual: tan# = y/x). It occurs|27|] that these data can be perfectly 



described by the solution of the classical equations of motion for the Abelian Higgs model 
(j78|). In Fig. |] by solid lines we show E z and kg for the Abrikosov string of the Abelian 
Higgs model ( |78| ) . The discretized version of the classical equations of motion is used in 



order to imitate the coarseness of the lattice used in numerical simulations [39] 



The parameters of the model corresponding to the best fit of the numerical data are: 
Q m /2n = 0.9519 ±0.0041, (83) 
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m v = 0.4522 ± 0.0206 = (1.0351 ± 0.0472) GeV , 
m H = 0.4747 ± 0.0600 = (1.0866 ± 0.1373) GeV . 



(84) 
(85) 



It appears that vector and scalar masses are equal to each other within the errors, thus 
we are at the so-called Bogomolny limit of the model. The Abrikosov string tension turns 
out to be 

y/a = 0.1808 ±0.0213 « (414.8 ± 48.9) MeV « 0.94^0^7^, (86) 

here asu(2) is the string tension of SU (2) gluodynamics. To get the dimensional quantities 
we have used y/asu{2) — 440 Mev. 



5.3 Anatomy of SU(2) Monopole 

At small values of the bare charge (at large values of /3) the compact electrodynamics 
is in the deconfmement phase, and gluodynamics is in the confinement phase. On the 
other hand, at large values of (3 in the maximal Abelian projection cos (pi (we use the 
parameterization ([40])) is close to unity (due to ( |38| ) and fl4"T|)), and the gluodynamic 
plaquette action \tr\Jp is close to the action of cQED, cosOp. Why the monopoles are not 
condensed in cQED and are condensed in gluodynamics if the actions of both theories are 
close to each other? It occurs that monopoles in gluodynamics have nontrivial structure 
and even near the continuum limit they differ much from the Abelian monopoles. The 
action of monopoles in gluodynamics is smaller than the action of monopoles in cQED. In 
Refs. g0[ it was found that the action of the non-diagonal gluons, , on the plaquettes 



near the monopole is negative, and the full non Abelian action, S su ^ = + S Abel , is 
smaller than the Abelian part of the action. Thus the action of monopoles in gluodynamics 
is smaller than that in Abelian theory and this is the possible explanation of the fact why 
the Abelian monopoles in gluodynamics are condensed at any value of the bare coupling. 

As we discussed in Section ^75] there exists the pure gauge field fl2"T|), (|28| ) which 
Abelian (A 3 /) component contains Dirac monopole. The lattice analogue of this field has 
zero non-Abelian action since the action of the Dirac string is zero on the lattice. And 



this configuration presumably corresponds [13] to the fields inside Abelian monopoles in 



gluodynamics. At large distances due to quantum fluctuations the fields around monopole 
may have a finite non-Abelian action density. 

Below we present the recent results |4l[] of numerical calculations in lattice SU(2) 



gauge theory which confirm the above discussion (see that paper for details). In Fig. [5| we 
show the dependence of S su ( 2 ' =< — S su ^ > on the half of the lattice spacing, a/2, 

in ferrmF]. The explanation of the scale of the horizontal axis is: < S^J^f 1 > is measured 
on the plaquettes which are faces of the cube dual to monopole, thus in a sense we are 
measuring the average field strength on the distance a/2 from the monopole center. 

The circles on Fig. |5] correspond to the calculation which take into account all monopoles, 
the squares correspond to monopoles taken from the percolating cluster (IR cluster). The 

2 We find the correspondence between the bare charge and lattice spacing fixing the value of the string 
tension a = 440 MeV and using the numerical data 142) for the string tension in lattice units, a ■ a 2 . 
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Figure 5: The dependence of S su ^ on a/2, the calculations are performed on the lattices 
up to 30 4 . 

results of the analogous calculation of the Abelian action near the monopole, S Abel =< S^- 
gAbei ^ are p resen ted in Fig. |6|. 

The results which follow from Figs. [|, ^| can be summarized as follows: 

• S su ^ at the center of the monopole belonging to the IR cluster is compatible with 
zero (at least it decreases when we approach the center). Thus at the center of 
monopole the fields are close to our prediction, eq. (|28|). The small value of S su ^ 
is the reason why IR monopole clusters percolate. Note that the percolation reflects 
the existence of the monopole condensate. 

• S 1 for monopoles belonging to IR and UV clusters is approximately the same. 
Thus there is no difference in the Abelian part of the monopole fields in IR and UV 
clusters. 

• S Abel increases when we approach the center of the monopole. The closer we are 
to the center (the larger fj) the larger lattice we have to use since the properties 
of the IR monopole cluster are strongly affected by the finite volume effects | ]32|j . 
This technical difficulty prevents us to decide whether the Abelian action density 
diverges or finite at the center of the monopole. 

Note that these facts (and also Figs. |5|, |^) are given for the action densities in lattice 
units, if some action density is constant in lattice units, it corresponds to 1/a 4 (Coulombic) 
behavior of this density in the continuum limit. 
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• all monopoles 

■ monopoles from IR cluster 
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Figure 6: The same as in Fig. || for S Abel . 
5.4 Conclusions ^ 4 

The results of numerical experiments in SU(2) lattice gauge theory confirm the validity of 
the dual superconductor model of the gluodynamic vacuum. The monopoles are condensed 
in the confinement phase, and already at the classical level are responsible for 94% of the 
SU(2) string tension. The structure of monopoles in gluodynamics is nontrivial, the 
non-diagonal gluons reduce the monopole action. 
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